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Abstract 

We construct explicit multivortex solutions for the first and second complex 
sine-Gordon equations. The constructed solutions are expressible in terms of 
the modified Bessel and rational functions, respectively. The vorticity-raising 
and lowering Backlund transformations are interpreted as the Schlesinger 
transformations of the fifth Painleve equation. 
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Motivation. Recently there has been an upsurge of interest in the complex sine-Gordon 
equation. Originally derived in the reduction of the 0(4) nonlinear a-model |I| and a 
theory of dual strings interacting through a scalar field [Q], this equation reappeared in a 
number of field-theoretic || and fluid dynamical 0] contexts. The equation was shown to 
be completely integrable , and the multisoliton solutions were constructed in a variety 

of forms, both over vanishing [|?],[8| and nonvanishing backgrounds PJlO|. The study of its 
quantized version started in [^,|TT],|l2| and received a new impetus recently |13] when it was 
realized that the complex sine-Gordon theory may be reformulated in terms of the gauged 
Wess-Zumino-Witten action and interpreted as an integrably deformed SU(2)/U(l)-coset 
model . 



The complex sine-Gordon theory can be conveniently defined by its action functional, 

d 2 x 
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The subscript 1 serves to distinguish this model from another integrable complexification of 
the sine-Gordon theory, the so-called complex sine-Gordon-2: 
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The latter system was derived in ref. [[TJ|] as the bosonic limit of a generalized supersymmetric 
sine-Gordon equation and, independently, in ref. fl6|]. Quantum mechanically, the above two 
complex sine-Gordon models were shown to be the only 0(2)-symmetric theories whose S- 
matrix is factorizable at the tree level [12[ . 



In all previous analyses the complex sine-Gordon equations were considered in the (1+1)- 
dimensional Minkowski space-time. In the present Letter we study these two models in the 2- 
dimensional Euclidean space. One reason for this is that they define integrable perturbations 
of Euclidean conformal field theories; more precisely, eqs.([L])-([|) arise as reductions of the 
SU(2) N gauged Wess-Zumino-Witten model perturbed by a multiplet of primary fields (by 
(fr^ and & 2 \ respectively) [nylT]]. They are closely related to important two-dimensional 
lattice systems, viz. parafermion models perturbed by the first and second thermal 
operators, respectively fl8| . 



Another motivation for studying solutions of the Euclidean complex sine-Gordon equa- 
tions comes from a remarkable similarity between eqs.(|I])-(@) an d several phenomenological 
lagrangians of condensed matter physics, in particular the Ginsburg-Landau expansion of 
the free energy in the theory of phase transitions, 
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and the energy of the Heisenberg ferromagnet with easy-plane anisotropy p0[ : 



E FM = J [(Va) 2 + sin 2 a(V/3) 2 + cos 2 a] d 2 x. 



(4) 



(Hence we will be using the words "action" and "energy" interchangeably in what follows.) 
To see that (0)-(@) are relatives of (f|), one writes ip = smae l/3 and ip = \/2 sin(a/2)e i/3 , 
transforming eqs.(|I]) and ([|) into 
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and 
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respectively. 

The Ginsburg-Landau free energy @ is minimized by the Gross-Pitaevski vortices orig- 
inally discovered in the context of superfluidity []19 |. These are topological solitons of the 
form ip(x,y) = $ n (r)e me , where $ n — > 1 as r — > oo. Although these important solutions 
were obtained numerically and in various asymptotic regimes, no analytic expressions for the 
Gross-Pitaevski vortices are available. Similarly, eq.(f|) is minimized by magnetic vortices 
PHI , and again, these are available only numerically. The aim of this note is to demonstrate 
that the Euclidean complex sine-Gordon equations also exhibit topological soliton solutions. 
Unlike the Gross-Pitaevski vortices and unlike their magnetic counterparts, the vortices of 
eqs.([TJ) and (0) can be found exactly, and in a closed analytic form. Consequently, the sig- 
nificance of the complex sine-Gordon equations on the plane stems from the fact that they 



provide a laboratory for studying analytic properties of vortices and their phenomenology 
in a wide class of condensed matter models. 

We construct these solutions in two different ways: (i) by means of an auto-Backlund 
transformation resulting from the spinor representation of the complex sine-Gordon theory, 
and (ii) via the Schlesinger transformation of the fifth Painleve equation, 
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which arises in a self-similar reduction of eqs.(|l|) and @. 

Vortices via Backlund transformation. The complex sine-Gordon-1 equation 
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admits an equivalent representation in terms of the Euclidean spinor field, $ = (u, v ) [JlO 



idu + v — \u\ 2 v = 0, 
idv + u — \v\ 2 u = 0. 



(9a) 
(9b) 



[Here d = d/dz, d = d/dz and z = (x + iy)/2.\ This is nothing but the Euclidean version 
of the massive Thirring model; the corresponding action functional has the form 
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Since as one can easily check both u and v satisfy eq.(g), the Thirring model (^) can be 
regarded as a Backlund transformation between two different solutions of eq.(§). Here we 
confine ourselves to multivortex solutions of the form if) = $ n (r)e me , where (r, 8) are polar 
coordinates on the plane and $ n (r) is a real function satisfying 
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Eqs.f]^) with u and v of the form u = — iQ n -\e t(jl ^ and v = $ n e me furnish an equivalent 
representation for eq. ([TT|) : 



- ^ + = (1 - (12a) 

^ + -$n = (1 - & n )*^l, (12b) 

dr r 

where <& n and satisfy eg. ([11]) with n and n' — n — 1, respectively. When n — 1, eq.( |12a| ) 
is solved by $ = 1 and eg. ( |12b| ) becomes a Riccati equation: 

$i + r" 1 $i = 1 - (13) 

This equation can be linearized by writing $! = S'/S, where S(r) satisfies the modified 
Bessel's equation of zero order: S" + S'/r — S = 0. Selecting S = Io( r ) gives the explicit 
form of the n — 1 vortex solution of the complex sine-Gordon theory: 

(14) 

Here /o( r ) an d A = -^o( r ) are the modified Bessel functions of zero and first order, respec- 
tively The vortex is plotted in Fig.l. 

With the solution $i at hand, eqs.([L2|) yield a recursion relation allowing us to construct 
solutions with vorticity n > 1 in a purely algebraic way: 
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In particular, the first two higher-order vortices (shown in Fig.l) are given by 
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where we have eliminated derivatives by means of the well known relation between the 
modified Bessel functions of different order: I n+ i + I n -\ = 2I' n . The asymptotic behaviour 
of the vortex with vorticity n is readily found from eq.( |12b| ): 

$n l —J n ~Zo7~ rJ n+2 + 0(r n+4 ) as r^O, (16) 
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$ n ~ 1 - — - — + 0(r~ 3 ) as r -> oo. (17) 
2r 8r^ 



One consequence of eq.([T7|) is that the energy of the vortices diverges [cf. eq.(^) below], 
similarly to the energy of the Gross-Pitaevski and easy-plane ferromagnetic vortices |T9| , |20|| . 
(Physically, this fact simply indicates that there is a cut-off radius in the system, for example 
the radius of the cylindrical superfluid container, or the distance between two adjacent vortex 
lines.) 

Bogomol'nyi bound. An important question is whether the vortex renders the action a 
minimum. Let n = 1 and rewrite eq.([l]) as 

Esg-i = J |00 + M 2 - 1| 2 Y^- 2 + J V ■ Ad\ (18) 

where A is a real vector field with components 

At = ln(l - |Vf) eySfArg^ + 2^; i = 1, 2, 

and ip = ipx + iip2- Assume our fields are such that \ip\ 2 < 1; then the first term in fll8|) 
attains its minimum at solutions to the "Bogomol'nyi equation" dtp = 1 — \ip\ 2 - This is 
exactly our eq. (prj|) with v = ip and u = —i; its vortex solution is given by eq.(^). The 
second integral in (|I8"D represents the divergent part of the action; it can be written as a flux 
through a circle of the radius R — > oo. Perturbing the vortex ip = &i(r)e 10 by a function Sip 
decaying faster than 1/r at infinity will not affect this part; the flux is uniquely determined 
by the vortex asymptotes: 

I A-ndl = 2n(2R-\nR-l) + 0(R~ 1 ). (19) 
Jc R 

Consequently, the n = 1 vortex saturates the minimum of the action in the class of functions 
with |V>| 2 < 1. 

The importance of the last inequality should be specially emphasized. Without the 
condition \ip\ 2 < 1 being imposed, one could construct a perturbation ip(x,y) of the vortex 
satisfying \ip\ = 1, = on some closed curve on the (x, y) -plane which does not enclose 
the origin. Taking then \ip\ 3> 1 in the interior of this contour, the action (|I]) could be made 
arbitrarily negative. 



It is interesting to note that the first-order equations (|) with generic u and v can also 
be interpreted as the Bogomol'nyi limit for some more general system with twice as many 
degrees of freedom. The corresponding action functional is 



E[u, v] = J f ^- ^ + 1 - \u\ 2 j d 2 x + j f ^- ^ + 1 - \v\ 2 j d 2 x + E Th} (20) 

where Etk is the Thirring action (0). Clearly, any solution to eq.@ is automatically a 
solution to the second-order system (|2"0|) . The action (p0|) can be written as 
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where Ai = ln(l — \v\ 2 )eijdj Aigv — ln(l — \u\ 2 )eijdjAigu. Assuming, again, that |m| 2 , |t>| 2 < 1, 
the lower bound of the action ( |20"| - [21| ) is saturated by solutions to eqs.([|). 

Some properties of the complex sine-Gordon vortices receive a natural interpretation 
when the equation is reformulated as a cr-model on a two-dimensional surface £ embedded 
in a three-dimensional space (n\, n 2 , n 3 ). The metric on E is ds 2 = da 2 + tan 2 ad/3 2 [see 
eq.(||)]. In order for £ to be smooth, the space (nj, n 2 , n 3 ) has to be pseudoeuclidean and 
the surface noncompact; in fact it looks like an asymptotically conical infinite bowl: 

ni + in 2 = tana e 4/3 ; 

1 i cos a 

riz = tanh q, q 



q ' ' (1 + cos 2 a) 1 / 2 

Here < a < 7r/2, < (3 < 2tc. In terms of n,, the lagrangian ([!]) reads 



E SG -i = j [(Vn 1 ) 2 + (Vn 2 ) 2 -(Vn 3 ) 2 + (l + n 2 + n 2 )- 1 ] d 2 x. 

As r — > oo, all three components of the vortex field, ni, n 2 and n3, tend to infinity. Conse- 
quently, the vortices map a noncompactified (x, w)-plane onto a noncompact surface — this 
accounts for their infinite energy. We also acknowledge the role of the condition \ip\ 2 < 1, 
which characterizes solutions admitting the a-model interpretation. 
Reductions to the Painleve-V. The transformation 

1 + W 
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reduces eq.([Ll|) to the fifth Painleve equation ([?]) with coefficients 

a = n 2 /8, (3 = -n 2 /8, 7 = 0, 5 = -2. (22) 
For 7 = c(l — a — b), where a 2 = 2a, b 2 = —2/3, and c 2 = —25, eq.(|7|) admits a reduction 



21[ to a Riccati equation 

W r = r-\W -\){aW + b) + cW. (23) 

The above relation between the coefficients is in place for n = 1; in terms of the vortex 
modulus $1, eg. (|23|) turns out to be nothing but our eq . (|13|) . Next, the Schlesinger trans- 



formations of the Painleve- V to itself [22,23] have the form 
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W r = r" 1 (W - l)(aW + b) + cW , , (24a) 

... 1 + W 
-W r = r- 1 (W-l)(aW + b)+cW —— . (24b) 

1 — W 

Here W and W satisfy eq.(0) with the coefficients (a,/3,7,5) and (a, (3,y, S), respectively, 
where a 2 = 2a, b 2 = — 2(3, 7 = c(b — a) , and 2d = a + b— 1 — 7/c, 26 = a + b — l + 7/c. With 
a, (3 and 7 as in eq . (|22|) , eqs. (|24]) amount to the vorticity-raising transformations (|T2|) . 

We conclude the discussion of the complex sine-Gordon-1 equation by mentioning that it 
would be natural to expect its vortex solutions (confined to a finite region on the plane) to 
arise as degenerate cases of its iV-soliton solutions |1(J (which have the form of N intersecting 
infinite folds). This kind of correspondence between two-dimensionally localized "lumps" 
and one-dimensional multisolitons exists, for example, in the Kadomtsev-Petviashvili equa- 
tion Surprisingly, the only two-dimensionally localized bounded solution resulting from 
the "degeneration" of the generic two-soliton solution of eq. (§) is discontinuous at the origin: 
ip = (X 2 — sinh 2 y)(X 2 + sinh 2 y) _1 . Here X + iY = e ta {x + iy), and a is an arbitrary 
constant angle. 

Vortices of the complex sine- Gordon- 2. The complex sine-Gordon-2 results from the 
variation of eq.(Q): 

V 2 ^ + ^| + - H 2 )(2 - H 2 ) = 0. (25) 



The multivortex Ansatz ip = $ n (r)e me = Ql/ 2 (r)e m6 takes it to 
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where a„ = and &„ 



-2n. Next, the substitution 



Q n = 2(l-W)- 1 

transforms eq. (|26|) to the Painleve equation (|7]) with coefficients 

a = 0, (3 = -2n 2 , 7 = 0, 5 = 2. 

This time, in order to construct the multivortex solutions we apply the Schlesinger trans- 
formation fl2~4f) twice. This leads to a recurrent relation 
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where and satisfy eq.(|26|) with the parameters (a k , b k ) and (a k -i,b k -i), respec- 

tively. Here a k _i = a k — 1 and b k _ 1 = b k — l. Starting with a trivial solution Q^> = 1 arising 
for a = —60 = n , an d using eq. (|27j) n times, we end up with a solution Q n = Q^~ n "> which 
satisfies eq.(|26|) with a n = and b n = —2n and the boundary condition Q n — * 1 as r — > 00. 
These solutions are given by rational functions; in particular, the first three multivortices 
(see Fig.l) read 



r 2 + 4' 



r 4 (r 2 + 24) ! 



^ 2 r s + 64r 6 + H52r 4 + 9216r 2 + 36864 ' 
Q 3 = r 6 (r 6 + 144r 4 + 5760r 2 + 92160) 2 ZV\ 

D 3 = r 18 + 324r 16 + 41472r 14 + 2820096r 12 + 114130944r 10 + 291962880O 8 
+ 50960793600r 6 + 611529523200r 4 + 4892236185600r 2 + 19568944742400. 



The energy of the complex sine-Gordon-2 vortices is logarithmically divergent. 



Concluding remarks. The Ginsburg-Landau expansion (|3|) is regarded as a central postu- 



this is different from fl3|) only in the third term which is small both when ip ~ and when 
l^l ~ ljVV" ~ 0. More importantly, the complex sine-Gordon models provide a unique 
opportunity for studying a number of analytic properties which are common to a wide class 
of vortex-bearing systems. These include the correct Ansatz for two spatially separated 
vortices, the vortex-phonon scattering matrix and so on; our present construction of coaxial 
multivortices is hopefully but a first step in this direction. Finally, one may see the complex 
sine-Gordon vortices as a starting point in the perturbative construction of the corresponding 
solutions of the Ginsburg-Landau and ferromagnet models. 

We are grateful to M. Bogdan, B. Dubrovin, B. Ivanov, A. Kapaev, P. Winternitz and 
the referee for useful remarks, and to N. Alexeeva and A. Harin for their numerical assistance 
at various stages of this work. One of the authors (I.B.) thanks S. Randjbar-Daemi for the 
hospitality at ICTP. This research was supported by the FRD of South Africa and URC of 



late in the phenomenological theory of phase transitions; however, for some systems eqs. ([[])- 
(0) may happen to provide a more adequate description. In fact, the difference is not as big 
as one might think. Assuming, for instance, \ip\ 2 < 1, eq.(|2]) can be rewritten as 




UCT. 



10 



REFERENCES 



[1] K. Pohlmeyer, Commun. Math. Phys. 46 (1976) 207. For the generalization of these 
results see J. Tafel, J. Math. Phys. 34 (1993) 1892, where the complex sine-Gordon 
equation is derived as a symmetry reduction of the self-dual Yang-Mills equations on a 
self-dual manifold. 

[2] F.Lund and T. Regge, Phys. Rev. D14 (1976) 1524 

[3] A.Neveu and N. Papanicolaou, Commun. Math. Phys. 58 (1978) 31 

[4] Y. Fukumoto and M. Miyajima, J. Phys. A: Math. Gen. 29 (1996) 8025 
Kulish, 

[5] F. Lund, Phys. Rev. Lett. 38 (1977) 1175; Ann. Phys. 115 (1978) 251 

[6] B.S. Getmanov, Pis'ma Zh. Eksp. Teor. Fiz. 25 (1977) 132 [JETP Lett. 25 (1977) 119] 

[7] H.J. de Vega and J.M. Maillet, Phys. Rev. D28 (1983) 1441 

[8] I.V. Barashenkov and B.S. Getmanov, Commun. Math. Phys. 112 (1987) 423 

[9] B.S. Getmanov, Teor. Mat. Fiz. 38 (1979) 186 [Theor. Math. Phys. 38 (1979) 124] 

[10] I.V. Barashenkov, B.S. Getmanov, and V.E. Kovtun, J. Math. Phys. 34 (1993) 3039; 
I.V. Barashenkov and B.S. Getmanov, ibid. 3054 

[11] H.J. de Vega and J.M. Maillet, Phys. Lett. B101 (1981) 302; J.M. Maillet, Phys. Rev. 
D26 (1982) 2755 

[12] G. Bonneau, Phys. Lett. B133 (1983) 341; G. Bonneau and F. Delduc, Nucl. Phys. B250 
(1985) 561 

[13] N. Dorey and T.J. Hollowood, Nucl. Phys. B440 (1995) 215; Q-H. Park and H.J. Shin, 
Phys. Lett. B 359 (1995) 125; T.J. Hollowood, J.L. Miramontes and Q-H. Park, Nucl. 
Phys. B445 (1995) 451 

11 



[14] I. Bakas, Int. J. Mod. Phys. A9 (1994) 3443; Q-H. Park, Phys. Lett. B328 (1994) 329 

[15] S. Sciuto, Phys. Lett. B 90 (1980) 75 

[16] B.S. Getmanov, Teor. Mat. Fiz. 48 (1981) 13 

[17] V.A. Brazhnikov, Nucl. Phys. B 501 [FS] (1997) 685 

[18] V.A. Fateev, Int. J. Mod. Phys. A 6 (1991) 2109 

[19] V.L. Ginsburg and L.P Pitaevskii, Sov. Phys. JETP 34 (1958) 858; E.M. Lifshitz and 
L.P. Pitaevskii, Statistical Physics, Part 2 (Pergamon Press, Oxford, 1980) 

[20] A.M. Kosevich, V.P Voronov and I.V. Manzhos, Sov. Phys. JETP 57 (1983) 86; V.G. 
Bar'yakhtar and B.A. Ivanov, Sov. Sci. Rev. A 16 (1992) 1; B.A. Ivanov and A.K. 
Kolezhuk, Low. Temp. Phys. 22 (1995) 275 

[21] N.A. Lukashevich, Diff. Eqns. 4 (1968) 732 

[22] B.I. Gromak, Diff. Eqns. 12 (1976) 519 

[23] A.S. Fokas, U. Mugan, and M.J. Ablowitz, Physica D 30 (1988) 247 

[24] L.A. Bordag, A.R. Its, V.B. Matveev, S.V. Manakov, and V.E. Zakharov, Phys. Lett. 
A63 (1979) 205 



12 



Figure Caption 

Fig.l The vortex solutions with n = 1, 2 and 3. 
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